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Abstract. The gravitational potential responsible for the lensing effect in SBS 1520+530 is studied over length 
scales from a few arc-seconds to a few arc-minutes. For this purpose, we use sharply deconvolved Hubble Space 
Telescope images in the optical and near-IR, in combination with ground based optical data obtained over a wider 
field-of-view. In particular, we have carried out a multi-color analysis in order to identify groups or clusters of 
galaxies along the line of sight. Photometric redshifts are measured for 139 galaxies unveiling significant excesses 
of galaxies 1.0' NW and 1.7' SW of the main lensing galaxy. The photometric redshift inferred both for the main 
lensing galaxy and for the galaxy concentrations is z=0.9^q'25- This is in rough agreement with the measured 
spectroscopic redshift of the main lensing galaxy, z=0.71 (Burud et al. 2002), suggesting that it is part of a larger 
group or cluster. We investigate the impact of including the galaxy cluster, first on the modelling of the lensing 
system, and second on the expected time-delay between the two quasar images. 

Key words, galaxies - clusters - gravitational lensing - multiple quasar system - individual SBS 1520+530 
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1. Introduction 

The study of multiply imaged quasars is one of the most 
promising way to measure cosmological parameters such 
as the Hubble parameter H (Refsdal 1964, Blandford & 
Narayan 1992). Indeed, H is related to two observables: 
the time delay between the light curves of the quasar im- 
ages and the mass distribution in the lens. Until recently, 
simple mass distributions were used to model multiply 
imaged quasars. However, the improvement of observing 
techniques has led to the discovery that most image con- 
figurations require more complex models, involving, for 
example, an external shear (Keeton, Kochanek & Seljak 
1997). This small but significant external perturbation to 
the main lensing potential often corresponds to the pres- 
ence of groups or even clusters of galaxies along the line 
of sight (see for example Keeton & Kochanek 1997, Burud 
et al. 1998, Morgan et al. 2001, Fassnacht & Lubin 2001). 
These additional structures must be taken into account in 
order to properly model the mass distribution and accu- 
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rately convert the observed time-delay into Ho. 
SBS 1520+530 is a doubly imaged quasar at z=1.855, with 
an angular separation of 1.568" (ct2ooo=15h 21m 44.83s , 
^20oo=+52° 54' 48.6"). It was discovered by Chavushyan 
et al. (1997) in the course of the Second Byurakan Survey 
(Markarian & Stepanian, 1983). Crampton et al. (1998) 
identified the lensing galaxy from near-IR adaptive op- 
tics images and proposed a simple lens model including 
exclusively the main lensing galaxy. They concluded how- 
ever that an additional external shear was mandatory to 
model properly the system. 

Deep ground-based and HST data are used in the present 
paper to map the mass distribution along the line of sight 
to SBS 1520+530. We then investigate the source of ex- 
ternal shear and propose a multi-components model for 
the total lensing potential, where the main lensing galaxy 
is a member of a larger group or cluster of galaxies. 



2. Observations - data reduction 

2.1. Ground-based observations 

The ground-based optical data consist of two sets of im- 
ages obtained during two observing runs, with the 2.56m 
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SBS 1520+53 (NICMOS2. F160W) 
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Fig. 1. Simultaneous deconvolution of the four HST/NICMOS2 frames obtained of SBS 1520+530. The left panel 
shows a standard combination of the data. The right panel is the result of the deconvolution, with a final resolution 
of 0.075". The pixel size on this 19"wide image is 0.0375". North is to the top, East is to the left. The labels follow 
the description in section 3.2. 



Nordic Optical Telescope (NOT) on the island of La 
Palma, Spain. 

SBS 1520+530 was observed with the High Resolution 
Adaptive Camera (HiRAC) on the night of February 25, 
1998 under photometric conditions and average seeing be- 
tween 0.8" and 1.2". The total exposure times were 90 
min in /, 40 min in R and 30 min in V. With a pixel size 
of 0.11", the useful field-of-view of the 2102 x 2052 de- 
tector is about 3.5' on a side. The images were reduced 
using sky flat-fields. Strong fringes (10%) are affecting the 
/-band images, but many dithered short exposures were 
combined into a high signal-to-noise map of the fringes 
which was then subtracted from the data. This procedure 
leads to final images with sky subtraction accurate to the 
percent level. 

Additional images were obtained on the night of 
February 16, 1998, with the Andalusia Faint Object 
Spectrograph and Camera (A1FOSC). The pixel size is 
0.188" on this instrument, which sets the field-of-view to 
6.5'. The observations were taken under 0.9" seeing and 
through thin cirrus. Total exposure times were 25 min in 
R and 30 min in V. These images were combined with the 
HiRAC (flux calibrated) data to obtain a deeper dataset 
but were also used to look for overdensities of galaxies over 
a wider field-of-view than that of the HiRAC data. 



2.2. Hubble Space Telescope observations 

Archival data of SBS 1520+530 were used to derive the 
precise astrometry of the system in the immediate vicinity 
of the lens. These data are part of a much larger public 
survey of gravitational lenses (PI: E. Falco) known as the 
Cf A- Arizona Space Telescope LEns Survey (CASTLES). 
NICMOS2 iJ-band (F160W) data were taken on July 20, 
1998, with a total exposure time of 47 min. At the date 
of observations, the NICMOS2 pixel size was 0.0759" x 
0.0753". 

Optical data were obtained with WFPC2 through the 
F814W filter (August 22, 1999), with a total exposure 
time of 27 min and through the F555W filter (November 
11, 1999), for a total exposure time of 35 min. 

3. Image deconvolution 

3.1. Image processing 

Several dithered images are available for each filter: 3 in 
F555W, 4 in F814W and 4 in F160W. These images were 
combined after cosmic rays and bad pixels rejection using 
standard IRAF procedures and were subsequently decon- 
volved using the "MCS" deconvolution algorithm (Magain 
et al. 1998). The algorithm can either deconvolve one sin- 
gle frame or deconvolve simultaneously several frames on 
the same object, combining them at the same time into 
one deep sharp frame. The algorithm also decomposes the 
data into a sum of point sources (here the quasar im- 
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Table 1. Information about the doubly imaged quasar and the lensing galaxies. The lines successively provide: (1) date 
of the observations, (2) magnitudes of the A and B images of the quasar (see designation in Fig. ||), coordinates relative 
to A, and spectroscopic redshift (Chavushyan et al. 1997), (3) flux ratio of the two quasar images, (4) magnitudes, 
coordinates relative to A and photometric redshift for galaxy L (with error bars at lcr), (5) similar data for galaxy M. 



ages and nearby stars) and extended sources (the lensing 
galaxy). One can therefore derive the photometry of the 
lens, decontaminated from the quasar light. 

We deconvolved the data in two ways. The F555W and 
F814W data were deconvolved only after combining the 
individual images, because of poor signal-to-noise ratio. 
For the more numerous and deeper F160W images, the 
individual data frames were simultaneously deconvolved in 
order to produce a high quality image of the system. This 
process is the same as described in Courbin et al. (1998) 
with a slight deviation from the standard application of 
this technique. In the standard application, each frame 
is seen as the convolution of a unique ideal frame with 
the PSF associated to each image in the data set. In the 
present case, the PSF available in the data (stars "SE" 
and "NW" in Fig. 1) are faint compared to the quasar 
images and we chose to compute one single PSF for all 
frames. This PSF is constructed from all stars available in 
the data set, i.e. two stars per frame, times four frames. 
This has the advantage of providing a high signal-to-noise 
PSF with very good rejection of cosmic rays and bad pix- 
els but it also implicitly assumes the temporal stability of 
the PSF over the period of observation. The quality of the 
residuals after deconvolution (see Courbin et al. 1998 for 
more details) remains good, indicative that our assump- 
tion is actually correct. 

3.2. Results 

The result of the simultaneous deconvolution of the 
F160W images is displayed in Fig. [|. Objects are labelled 
as in Crampton et al. (1998). The two PSF stars used for 
the deconvolution are the ones labelled "SE" and "NW" . 
Note the elongated lensing galaxy L appearing between 
the two quasar images. Two other objects are also visible 
in the immediate vicinity of the lensed source (labelled 
galaxies M and N in Fig. |l]). Although at the confidence 
limit after the deconvolution process, the extended plumes 
to the NW of the lensing galaxy L might trace a prior in- 
teraction with galaxy M. However we shall ignore this fea- 
ture in the subsequent modelling. Table [l] summarizes the 
main properties of the quasar images and of the lensing 



galaxy. As the lensing galaxy appears to be essentially el- 
liptical, models are particularly sensitive to its orientation 
and ellipticity: hence particular attention was paid to the 
precise determination of these parameters. We used the 
Im2shape software (Bridle et al. 2002) to determine accu- 
rate values of the PA and of the ellipticity of the galaxies L 
and M. We measure an ellipticity of 0.50±0.02 for galaxy 
L, the ellipticity being defined as e = (a 2 — b 2 )/(a 2 + b 2 ). 
Its major axis is oriented at PA= —23° ± 7°. Galaxy M 
is located at 2.56" to the NE of the quasar image A and 
is likely to play some role in the lens modelling. It is al- 
most circular (e=0.04±0.02, PA= -40° ± 10°). Galaxy N 
is too far away from the main lensing galaxy L (~7") to 
play an important role in the modelling, at least if realistic 
masses are considered. The coordinates of the quasar im- 
age B, of the lensing galaxy L and of galaxy M, relative to 
the quasar image A, are provided in Table [j]. Photometric 
redshifts have been derived for the different objects in the 
immediate vicinity of SBS 1520±530. These redshifts were 
estimated with Hyper^ (Bolzonella et al. 2000), as were 
tho se of the intervening galaxy cluster/group (see Section 
4.2 for complete details). Results from Hyper z consist in 
two different redshift estimates. The first redshift corre- 
sponds to the best fit (smallest x 2 ) obtained after testing 
a wide variety of galaxy templates at different redshifts. It 
is called the "photometric redshift" , z p hot- The second red- 
shift corresponds to a weighted mean redshift z wm , repre- 
sentative of the "most frequent" redshift allocated to the 
galaxy, while testing a family of galaxy templates over 
a range of redshifts. The latter redshift, z wm , is less af- 
fected by degeneracies due to redundant spectral features 
in galaxy templates. We obtain, for galaxy M, Zph t= .89 
-o 26 an d z «jm=0.57, and for the main lensing galaxy L, 
z p hot0 -88 1q °6 an( i z Jum=0.71. At the time we were per- 
forming this analysis, a spectroscopic redshift of galaxy L 
was measured, z=0.71 (Burud et al. 2002). These three 
values are consistent within the error bars. 



1 The original code and user's manual are available at 
http://webast.ast.obs-mip.fr/hyperz 
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Fig. 2. The A1FOSC i?-band image (5.7' x 6.0') around SBS 1520+530. The isodensity contours outline regions with 
constant number of galaxies per square arc-minute. Notice the concentrations of galaxies to the NW and SW of the 
lensed quasar (indicated by an arrow). 



Fig. 3. Similarly to Fig. ||, isodensity contours deter- 
mined for the 2.8' x 2.7' HiRAC field (I-band) around 
SBS 1520+530. The edges of the overdensity of galaxies 
already seen in the A1FOSC image, are also detected in 
the NW and SW corners of the frame. 



4. Detection of an intervening galaxy 
cluster/group 

4.1. Photometry of the field: an overdensity of galaxies 

The photometry of all galaxies detected in the field- 
of-view near SBS 1520+530 was performed using the 
SExtractor 1.2 software (Bertin & Arnouts 1996). No ob- 
vious defects were present on the detector chip. Yet the 
central part of the field (20" in diameter) was masked, as 
it is saturated by a very bright star. 

The presence of an intervening galaxy cluster/group 
along the line of sight to SBS 1520+530 was first inves- 
tigated on the wide field dataset provided by A1FOSC. 
In the i?-band, we have detected 620 galaxies over a 6.5' 
field-of-view after removal of the stellar objects. A stel- 
lar object was defined as an object brighter than R = 19, 
with a FWHM smaller than 1.0" (the seeing being 0.9" ) 
and with elongation (defined as a/b) smaller than 1.1. 
Isodensity contours were obtained from the A1FOSC im- 
age, for all galaxies in the range 22 < R < 25. This 
eliminates low redshift "normal" galaxies but may still 
include some contamination by faint dwarf galaxies which 
cannot be identified using only V, R and / photometry. 
Following these criteria we have displayed the i?-band im- 
age in Fig. ^| where a significant overdensity of galaxies is 
seen close to the centre of the field. The projected density 
of galaxies is a factor two larger in the densest parts than 
in the rest of the field. The shape of the overdensity is 
irregular and even double-peaked, with one peak located 
at 1.0' NW of SBS 1520+530 and a second peak at 1.7' 
to the SW. These overdensities also appear in the smaller 
ficld-of-view of the HiRAC /-band image, applying simi- 
lar selection criteria, i.e. rejecting the stellar objects and 
considering only the objects with 21 < / < 24 (see Fig. ||). 

4.2. Photometric redshift analysis 

Multi-color information is available in a 2.8' wide field 
around SBS 1520+530 which corresponds to the HiRAC 
ficld-of-view. This is sufficient to recognise the concentra- 
tion of galaxies detected in the A1FOSC image and to infer 
photometric redshifts, based on their V, R and / magni- 
tudes. A total number of 139 galaxies are detected in the 
three filters simultaneously The magnitudes and their er- 
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Fig. 4. Relative distribution of galaxies with two different 
redshift estimators, for the HiRAC field-of-view around 
SBS1520+530 (solid line). Top panel: distribution for the 
photometric redshifts (z p ^ ot ). Bottom panel: distribu- 
tion for the weighted mean redshifts (z wm )- In each panel 
we have also plotted the distribution expected for a refer- 
ence field of the same size and containing no cluster/group 
(dashed line). 



ror bars are used as the input to Hyperz (Bolzonella et 
al. 2000), in order to estimate redshifts. Hyperz makes a 
comparison between the photometric spectral energy dis- 
tribution (SED) of the observed galaxy and those obtained 
from a set of reference template spectra, using the same 
photometric system. Reddening is also taken into account: 
we used the Calzetti law (Calzetti et al. 2000). The lim- 
iting magnitudes, corresponding to ~80% completeness, 
are 25 in V, 24 in R, and 23 in /. This allows to detect 
objects up to z^l.2. The entire library of template spectra 
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Fig. 5. Composite image corresponding to the HiRAC field-of-view (2.8' on a side taking into account dithering). 
The V, R and I images have been used to build this "true-color" image. In this image, violet galaxies (large blue 
circles) are allocated redshifts in the range 0.1<z<0.3 by Hyperz while the red ones (small yellow circles) are in the 
higher redshift bin 0.7<z<1.0, as the main lensing galaxy. 



proposed in Hyperz has been considered, encompassing a 
variety of star formation histories. 

The most likely redshift is estimated for each galaxy 
through a \ 2 minimization. As the distribution of % 2 ver- 
sus redshift is not symmetrical around the best redshift 
value, the weighted mean redshift z wm can also be a rel- 
evant and interesting parameter. In Fig. ^, we have plot- 
ted the relative distribution of galaxies as a function of 
redshift, using both the photometric redshift z p hot (upper 
plot) and the weighted mean redshift z„ (lower plot). 
The distributions are shown along with that of a simu- 
lated field with the same mean density of galaxies, but no 
galaxy clustering. In the simulation, the redshift of galaxy 
formation is set at Zf=7, and the cosmological parameters 
are Ho=50 km.s _1 .Mpc _1 and Qq=1, Two redshift bins 
show an excess of galaxies, at z^0.2 and at z~0.9. These 
two concentrations in the redshift space appear in both 
types of distributions, whether z p hot or z wm is considered 
(see the two panels in Fig. ||). Moreover, the shapes of the 
distributions do not change if a redshift limit is imposed 
in Hyperz, for example at 1.0 or 1.5. Therefore, the galaxy 
excesses at z~0.2 and z~0.9 are not related to the choice 
of the redshift limit. 

The color image on Fig. || illustrates this result as well: 
the galaxy population with violet colors corresponds to 
galaxies around z=0.2 while the galaxy population with 
red color corresponds to galaxies around z=0.9. 

As the galaxies contributing to the red population 
(redshift peak around 0.9) are also those which build up 
the overdensities seen in Fig. | and |, we deduce that the 
cluster/group located to the West of the quasar images 
is indeed at redshift z=0.9lg'25 ficr error bars). We can- 
not exclude that there might be a real difference between 
the redshift of the main lensing galaxy (z=0.71) and the 
redshift of the cluster/group. However, we consider in the 
modelling (see section ||) that they are both at z=0.71, as 
this is still consistent with the error-bars. 

5. Lens modelling 

Analytical modelling of SBS 1520+530 was performed us- 
ing the Lenstool software developed by Kneib et al. (1993). 
This software computes the lens mass distribution neces- 
sary to reproduce the configuration and flux ratios of mul- 
tiply imaged quasars. In the present case, multiple com- 
ponent lenses can be used, each one being a truncated 
Pseudo-Isothermal-Elliptical-Mass-Distribution (PIEMD , 
Kassiola & Kovner 1993, Kneib et al. 1996, Hjorth & 
Kneib 2002). Given the observational constraints avail- 
able we explored a number of lens configurations, involv- 
ing a single lensing galaxy (model L), two galaxies (model 
L+M), and two galaxies plus a galaxy cluster (model 



L+M+C). Obviously, the galaxy cluster is represented by 
a circular PIEMD mass distribution. 

The constraints available are the positions of the 
quasar images and their flux ratio in the F160W band 
(Table |l|), the positions of galaxies L and M and their PA 
(9) and ellipticities (e) with error-bars given in section || 
and the position of the cluster. The error-bars on PA take 
into account that the misalignment between dark-matter 
halo and light of a galaxy could be of 10° (Kochanek 
2001), so that the galaxy mass distribution could have an 
orientation different from that traced by the light. 

Given the results and uncertainties from the photomet- 
ric redshift analysis (see section 3 and 4), a same redshift 
can be considered for all model components. We take this 
redshift to be that of the main lensing galaxy, z=0.71. 

5.1. Model L 

We first run a model considering only the lensing galaxy L. 
The fitted parameters, r cut (L) and ctl are given in Table § 
The quantities and 6l, as well as the position of the 
lens, are fixed to their observed values within their respec- 
tive error bars. In this case, neither the flux ratio nor the 
relative position of the quasar images are recovered, as 
shown by the large value of the \ 2 ■ This strengthens the 
suspicion that galaxy L is not the only intervening lens. 

5.2. Model L+M 

Therefore, in a second step we consider that galaxy M 
(the nearest galaxy to galaxy L) contributes to the total 
lensing potential as well. Again the quantities ex, €m, 8l, 
6m , and the positions of the two galaxies are fixed to their 
observed values within their respective error bars. In this 
case we recover the flux ratio, but not the geometrical 
configuration of the quasar images. 

5.3. Model L+M+C 

Therefore, additional mass is necessary to fit both the 
configuration and flux ratio of the quasar images. A 
contribution from a galaxy cluster is proposed and intro- 
duced in the model. We first center the cluster between 
the two peaks detected in the galaxy density map (see 
Fig. ||), but we fail to obtain an acceptable \ 2 n t to the 
data. The best x 2 values are obtained for a centering of 
the cluster component on the NW peak of the galaxy 
overdensity. Doing this, we obtain the results summarized 
in the last line of Table |2[ Adding the SW peak of the 
galaxy overdensity as a forth PIEMD still provides a good 
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Table 2. Summary of the models studied. The quantity e is the ellipticity of the object, 9 is its PA, a its velocity 
dispersion and r cut its cut radius. Column (1) indicates which components have been considered (galaxy L, galaxy M, 
cluster C) and the indice or 1 represents the number of degree of freedom of the model. Column (2) provides the 
derived values for the fitted parameters. Column (3) summarizes the fixed parameters: in all cases the positions of the 
components were fixed to their observed values; the ellipticities for galaxies L and M where fixed to their observed 
values; when fixed, 6l is -23°and 8m is -40°; r cut (L)= 51 kpc, r cut (M)=20 kpc and r cut (C)= 600 kpc. Column (4) 
provides the \ 2 value. Column (5) gives the time-delay prediction expressed in days, for Ho=65 km.s .Mpc f2=0.3 
and A=0.7. Finally, Column (6) displays the mass-to-light ratio value A4/£j computed within the Einstein radius for 
galaxy L (solar units, see section ||). 



fit. Indeed, being further away, its contribution is minor 
in the lens model and its unique impact is to reduce the 
time-delay by a few days. 



5.4. Inferences from the lens modelling 

The mass to light ratio Ai/£x of galaxy L computed 
within the Einstein radius (r^ ^156 kpc or 0.7") for 
the different models is given in Table ||. The luminosity 
Cx was determined directly from the observed F160W 
photometry, as the F160W filter for a galaxy at z ~ 0.9 
roughly correspond to the I-band filter at z=0, and thus 
no color correction was applied. For an elliptical galaxy 
C-b/ C-i ~l-9 therefore the derived M/Cb for the galaxy 
L, is consistent with the range of usual M/Cq ratios 
found in lens systems (Keeton et al 1998). 

The time-delay predictions given in Table || were 
computed for the various models of the lens, using 
H =65 km.s -1 . Mpc -1 Q=0.3 and A=0.7. Changing the 
cosmology to fi=l and A=0, with the same value for 
Ho, decreases the time-delay to 85 days for the L+M+C 
model. Changing the value of Ho from 65 km.s -1 . Mpc -1 
to 50 km.s -1 . Mpc -1 with Sl=0.3 and A=0.7, increases the 
time-delay from 100 days to 125 days for the L+M+C 
model. 

6. Concluding remarks 

We have studied the doubly imaged quasar SBS 1520+530 
using both HST and ground based data, in order to map 
the projected mass distribution of the lensing system at 
small and large scales. 

First, by deconvolving the HST/NICMOS data, we 
have accurately measured the shape and brightness of the 
main lensing galaxy. We detect a weak extension to the 
NW of the lensing galaxy (Fig. 1). This feature does not 
seem to be a residual of the deconvolution process. If it is 
a genuine feature, it may have some influence on the lens 



model, considering its proximity to image A of the quasar. 
However until it has been confirmed, we ignore it for the 
modelling. 

Second, we have mapped the relative surface density 
of faint galaxies around the quasar images and identified 
an overdensity of objects to the West. Photometric red- 
shifts were estimated for the main lensing galaxy L as 
well as for all the galaxies in the 7.5 (') 2 surrounding the 
lensed quasar. This analysis reveals a concentration in red- 
shift space around a mean value z=0.9^o25 anc > confirms 
the 3D-reality of the observed projected overdensity of 
objects: we consider it to be a galaxy cluster. As a con- 
sequence, it is probable that the main lensing galaxy L 
(measured spectroscopic redshift z=0.71) is a member of 
this cluster (measured photometric redshift z=0.9lo 35)- 
Following this argument and considering that the spec- 
troscopic redshift measurement is more accurate, we have 
built our lens model with a unique redshift value, z=0.71, 
for the main lensing galaxy, for galaxy M and for the clus- 
ter. 

Finally, the presence of a galaxy cluster to the W of 
SBS 1520+530 is also suggested by the lensing model 
analysis which does require such a contribution. If not 
taken into account, it is impossible to reproduce the ob- 
served configuration and flux ratio of the double quasar 
images. This argument is valid as long as light traces 
mass, which seems to be the case for most lensing galaxies 
(Kochanek 2001). 

An independent and complementary way of probing 
the line of sight mass distribution is to measure the time- 
delay induced by the lens system (also necessary to con- 
strain Ho). Therefore, it is interesting to predict the time- 
delay corresponding to the lens model derived from the 
lensed quasar image configuration. We have made time- 
delay predictions for the 3 different lens models analysed 
in this paper, to be compared with the observed time- 
delay, when it will be available. 

Yet, further detailed observations are needed to im- 
prove the understanding of the lensing system and en- 
able to use it as a cosmological probe. SBS 1520+530 
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will be observed by the X-ray satellite Chandra (PI: 
G.Garmire): this observation should confirm the existence 
of the cluster in a way similar to the quadruply lensed 
quasar RX J0911.4+0551 (Morgan et al. 2001). Deep 
Keck/Gemini multi-object spectroscopy would be invalu- 
able to settle the redshift of the detected cluster/group 
and provide an estimate of its mass. 
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